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Abstract 

Motivated by recent CLEO measurements of the B —> rf K decay, we evaluate 
gluon/charm content of the rf meson using the interacting instanton liquid model of 
the QCD vacuum. Our main result is (0| g 2 ’ f abc G^ v G\, a G c ai f\rl) = (2.3 + 3.3)GeV 2 x 
(0| g 2 G c f JU G c f LV \rf). It is very large due to the strong field of small-size instantons. 
We show that it provides quantitative explanations of the CLEO data on the B —> 
rf K decay rate (as well as inclusive process B —> rf + X), via a virtual Cabbibo- 
unsuppressed decay into cc pair which then becomes rf. If so, a significant charm 
component should be present in other hadrons also. In particular, we found a large 
contribution of the charmed quark in the polarised deep-inelastic scattering on a 
proton. 


1.Instantons of a small size (p ~ 1/3 fm) are known for long time to be a very 
important component of the QCD vacuum (T|. In general, their fields are respon¬ 
sible for a scale 1 GeV which restrict perturbative QCD from below, and effective 
hadronic Lagrangians from above. Because of fermionic zero modes, they play es¬ 
pecially important role for light (u,d,s) quark physics (for recent review see JH). It 
was nevertheless believed that they are irrelevant for charm-related physics: and 
indeed, the instanton-induced spin-dependent and independent potentials between 
heavy quarks are small compared to standard confining-plus-perturbative one. How¬ 
ever, as we show in this paper, the situation is reversed for virtual cc pairs: they 
can only appear due to the strongest gluonic fluctuations in vacuum, and those are 
instantons. (In fact, the gluonic fields in the centre of relevant instantons is so large, 
that one may even question whether gG/m 2 is a good expansion parameter.) 

The way to see this is to look at the charm component in hadrons with different 
quantum numbers. The object of this paper, rf , is long known to play a very special 
role in QCD: separated by a large gap from other pseudo-scalars (the Weinberg’s 
U(l) problem |3], Q) it serves as a screening mass for the topological charge (see 
recent detailed discussion in pj). Thus testing whether the high dimension gluonic 
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operator does or does not couple strongly to the rj we are actually testing whether 
the strongest vacuum fluctuations do or do not possess the topological charge. No 
effect of such magnitude should exist e.g. for vector mesons: and indeed, the em¬ 
pirical Zweig rule is very strict in vector channels, allowing only tiny flavor mixing. 

2. Recently, CLEO collaboration has reported || measurements of inclusive and 
exclusive production of the rf in B-decays : 

Br(B -*• rf + X ; 2.2 GeV < E r/ i < 2.7 GeV ) = (7.5 ± 1.5 ± 1.1) • 10 " 4 , (1) 

Br(B -> rf + K ) = (7.8±|;| ± 1.0) • lO " 5 . (2) 

Simple estimates 0 show that these data are in severe contradiction with the stan¬ 
dard mechanism, the 6 -quark decay into light quarks, because Cabbibo suppression 
factor V u b leads to numbers which are by two orders of magnitude smaller than the 
data (both the inclusive and exclusive cases). Alternative mechanism, suggested in 
0 is based on the Cabbibo favored 6 —> ccs process, followed by a transition of 
virtual cc into the rj. The latter transition may be possible, provided there exist 
large intrinsic charm component of the rj'. Its quantitative measure can be expressed 
through the matrix element 

(0|cvysc| rj'(q)) = ifjfqp . (3) 

and one needs /(-, « 140 MeV in order to explain the CLEO data, see 0. This 

value is surprisingly large, being only a few times smaller than the analogously 
normalised residue ( 0 |c 7 At 75 c| 77 c (g)) = if Vc qji with f Vc ~ 400 MeV known experi¬ 
mentally from the ij c —> 77 decay. 

3. Because the c-quark is heavy, it may only exist in the rj' in a virtual loop,and 
its contribution can be evaluated in terms of gluonic fields. Taking the divergence 
of the axial current in Eq.([3j) one gets 

& = -^2-(0|2m c ci7 5 c+ ^-G^G^rj') . (4) 

v m, ;/ 4-7T 

which can be further simplified by the Operator Product Expansion in inverse pow¬ 
ers of the c—quark mass 

2m c ci lb c = ~^G, U G, U - ^ ^ff abc G%G b ua G c ^ + 0{G^m\) (5) 

(see the appendix in 0 for a detailed derivation of this result. Further terms in 
expansion (j5) are neglected in what follows.) Thus the problem is reduced to the 
matrix element of a particular dimension -6 pseudo-scalar gluonic operator: 

4 ° = . ( 6 ) 

T] c 

4. The magnitude of the matrix element (p) was related 0 to the vacuum ex¬ 
pectation value of similar operators: 

Ac) _3_1 ( 9 3 G 3 )ym m 

v ' ^bmtiOl^G^G^W) ' 
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where (G 3 ) should be evaluated in pure gluodynamics, not QCD. Unfortunately, 
only indirect order-of-magnitude estimate for this latter quantity was given, and 
thus in 0 rather wide range of values was given = (50 -f- 180) MeV. 

5.We have performed direct calculation of this quantity using the Interacting 
Instanton Liquid Model (IILM). In its present form, this model takes into account 
instantons coupling to light quarks to all orders in t’Hooft effective interaction, 
which was shown to be crucial for rf physics. It has correctly reproduced multiple 
mesonic/baryonic/glueball correlation functions, and also has an increasing direct 
support from lattice studies of instantons (see HI). 

The calculation is based on numerical evaluation of the following two-point Eu¬ 
clidean correlation functions 

I< 22 (x) = (0\g 2 G; u Gl u (x)g 2 G^Gl u m0) ( 8 ) 

K 23 (x) = (0|<7 2 G^G^(x) 5 3 / abc G^Gt A G^(0)|0) (9) 


K 33 (X) = <0| g 3 f abc G; u Gl x G c Xll (x) g 3 f abc G; u G b uX G c Xli mo) (10) 

Studies of K 2 2 (x) has been made previously | 8 |, where it was demonstrated that in 
the “unquenched” ensemble of instantons with dynamical quarks the non-perturbative 
part change sign at distances x > 0.6/m, displaying a “Debye cloud” of compen¬ 
sating topological charge. It is identified with the rf contribution, and lead to an 
estimate 

<0| g 2 G%G%\ V ') = ^U'rn%, « 7GeD 3 (11) 


which agrees reasonably well with other estimates in literature. In this formula 
we have expressed matrix element (0 in terms of the standard parameter /„/ ~ 
85MeV which is defined as follows 


<°|4 msQiW) 

V 8 i=u,d,s 


Ifr]' Qfi- 


Using an anomaly in the chiral limit, m u = m^ = m s = 0 we arrive to ©■ 

6 . We have calculated the correlators mentioned by numerical simulation, using 
of the ensemble 16 instantons and 16 anti-instantons, put into a box 4x 2 3 /m 4 , with 
(without) dynamical quarks [p] Unfortunately, the propagation of the gluons in the 
background non-perturbative fields of instantons was not studied in such details as 
for light quarks, and so far we do not have the gluon propagator program which could 
be used for all distances. At small x purely perturbative results (e.g. K^^ix) = 
384<7 4 /7t 4 x 8 ) dominate, while the non-perturbative fields can be included via the 
operator product expansion (see e.g. 0 -i)- At large x we would argue below that 
(at least with dynamical quarks) the non-perturbative fields dominate. 

The quantity ° (|) can be obtained from the correlation functions (| 8 |, ||, 

0 ) 0 : 


f? vj m 2 c 

frf 


K 23 {x -> 00) _ / K 33 (x -> 00) 
K 22 (x —> 00) V K 22 (x —^ 00) 


( 12 ) 
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It is expected that at large distances the contribution to two other correlators would 
also be dominated by non-perturbative field of the instantons. If so, one has a simple 
estimate for the ratio of matrix elements 


(0\g 3 f abc G^Gl a G c atl W) 12, 1 


(0\g 2 G“ G% W 


= T <-> 

5 p z 


(14-1.5 )GeV 2 


(13) 


Two numbers given here correspond to averaging over instanton size distribution for 
two variants of the instanton-anti-instanton interaction, the so called “streamline” 
and “ratio-ansatz” ones, and indicate the sytematics involved. The latter (giving 
smaller average size and larger number above) should be considered preferable, 
because it better agrees with the size distribution directly obtained from lattice 
gauge field configurations, see discussion in (^]. (Recent measurements [12] using 
refined “inverse blocking” method has found somewhat smaller instantons than 
others, but those seem to belong to correlated instanton-anti-instanton pairs, which 
would not contribute to the compensating Debye cloud we look for.) 

In our measurements of K23, K33 both ratios entering (12) were found to stabilize 
at large enough x> 0.8/m at the same numerical value. We take it as an indication 
that vj contribution does in fact dominate, although we were not able to see that 
all correlators fall off with the right mass (T^j . Numerical values of the ratios about 
(1.5 4- 2.2 )GeV 2 , for two ensembles mentioned. The nurners are somewhat larger 
than in (|l^) because the second operator in the correlator makes it more biased 
toward smaller instantons. 

Proceeding to final result, we have to look at radiative corrections. The experi¬ 
mental number mentioned above is defined at the scale pi ~ mj?, which is different 
from that obtained in the instanton calculation. In the latter case the charge and 
fields are normalized at p?, ~ gG where G is the typical gauge field at the points 
which contribute the most to the correlators. Two scales are not too far apart nu¬ 
merically P 2 ~ (0.54- l)GeV 2 , but the anomalous dimension of the g^GGG operator 
114] is large, and it leads to correction 


f^\pi ~ m c ) = [ 


a, 




Oi s (p2) 


l-5/,y V 2 ) 


(14) 


Here we use m c (pi — m c ) ~ 1.25 GeV for the numerical estimates. This concludes 
our derivation of the parameter 


r( c ) 

-j— ~ 0.85 4- 1.22, (15) 

Jr/ 

where the second value is preferable, see above. We present our final result as 
the ratio instead of the absolute value of because most systematic 

errors are gone for the ratio. The final uncertainty in eq. © comes from the 
systematic errors of the instanton model, which can be judged from comparison of 
the instanton size distribution or the scalar glueball size to corresponding lattice 
results. It will certainly be soon reduced by on-going works. Finally, we compare 
it with the “experimental” value needed to explain CLEO measurements ([|), and 
conclude that our result obtained in the instanton liquid model agrees with it, inside 
the uncertainties. 
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7. The next logical question to ask is whether the unexpectedly large gluon/charm 
content of rf has profound consequences outside of B physics, for other hadrons. 

One point we want to make is that that it seems now likely that understanding of 
the spin problem of the nucleon cannot be done without its “intrinsic charm” as 
well[0]. Relevant matrix element 

(N\c'y fl 'y 5 c\N) = N'y^N (16) 

could be generated by the rf “cloud” inside the nucleon. Assuming now the rf dom¬ 
inance in this matrix element 0] one could get the following Goldberger-Treiman 

type relation |l5fl g'^ = 2 Ai v fJv'NNf^ ■ Although the precise value of g^NN is un¬ 
known, and phenomenological estimates of the coupling vary significantly g^NN = 
3-7@ , it leads to 

(N\c'y ll 'y 5 c\N) = (0.2 Q.5)N'y fi 'y 5 N (17) 

which is comparable to the light quark contribution! We plan to calculate and 
grfNN in the instanton model as well. Lattice determination of all those quantities 
would be more than welcome. Ultimately, the contribution of the charmed quarks 
in polarized deep-inelastic scattering may be tested experimentally, by tagging the 
charmed quark jets. 

8.It is by now widely known that the Zweig rule is badly broken in all scalar/pseudoscalar 
channels, and that (rather large) mass of the rf is in fact due to light-quark-gluon 
mixing. Furthermore, all these phenomena are attributed to instantons. In this 
work we have found that similar phenomena are even more profound for larger- 
dimension (multi-gluon) operators as well. Moreover, the flavor mixing includes 
also a significant fraction of cc in rf. Perhapse it is not so surprising qualitatively: 
but the fact that one can actually quantitatively calculate these matrix elements 
and quantitatively compare it to real data is still rather amasing. 
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